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Abstract 

We study the Gross-Pitaevskii functional for a rotating two-dimen- 
sional Bose gas in a trap. We prove that there is a breaking of the 
rotational symmetry in the ground state; more precisely, for any value 
of the angular velocity and for large enough values of the interaction 
strength, the ground state of the functional is not an eigenfunction 
of the angular momentum. This has interesting consequences on the 
Bose gas with spin; in particular, the ground state energy depends 
non-trivially on the number of spin components, and the different 
components do not have the same wave function. For the special case 
of a harmonic trap potential, we give explicit upper and lower bounds 
on the critical coupling constant for symmetry breaking. 



1 Introduction 

We consider the Gross-Pitaevskii (GP) theory of a rotating two-dimensional 
Bose gas in a trap. The Bose gas is described by a single function on R^, 
the wave function of the condensate. It is confined in some trap potential 

*Erwin Schrodinger Fellow. On leave from Institut fiir Theoretische Physik, Universitat 
Wien, Boltzmanngasse 5, 1090 Vienna, Austria 
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V, and rotates around the origin at an angular velocity fi. The strength of 
the interaction between the particles is measured by the positive parameter a 
appearing in the GP functional ( |1.3| ) below. It is related to the particle num- 
ber A^, the scattering length of the interaction potential and the average 
particle density p via 

a=- — ^. (1.1) 

Minimization of the GP functional is supposed to describe the physical prop- 
erties of rotating Bose gases at very low temperatures, as considered in recent 
experiments |]T], |^. These show various interesting properties, in particular, 
the appearance of multiple vortices and a resulting breaking of the rotational 
symmetry. There have been a lot of theoretical investigations on these phe- 
nomena (see, e.g., the review article 0), based on the GP approach, either 
using numerical methods or various simplifying approximations, but a proof 
that the Gross-Pitaevskii functional indeed captures all these features is still 
missing. 

We shall not be concerned here with the derivation of the GP func- 
tional from the basic quantum mechanical A^-particle Hamiltonian. For non- 
rotating systems, i.e., f2 = 0, this has been achieved in 0, |] (see also P). 
However, the methods used there allow no simple generalization of these 
results to the rotating case. 

We shall now describe the setting more precisely. We denote by (r, y?) 
polar coordinates for x = {x,y) G M^. In these coordinates, the angular 
momentum is given by L = —id /dip. For 

Ho = -A-nL + V{r) (1.2) 

and (j) G Q{Hq) fl L^(R^,(i^x) define the Gross-Pitaevskii energy functional 
by 



= (0|ifo0) + a J |0(x)|Vx. (1.3) 

(Here Q{Ho) denotes the quadratic form domain of Hq, and (-I-) denotes the 
standard inner product on L^(R^)). The parameter a is non-negative, and 
without loss of generality also i7 > 0. We assume that V G Lj^^(]R^) is a 
positive radial function with the property that 

\/(r) > f]V/4 - (1.4) 



R. SEIRINGER, ROTATING BOSE GAS 



3 



for all fl in some non-zero interval [0,fic), with < oo. We take [0, ^c) 
to be the largest such interval, allowing it to be the whole half line; i.e., 
we allow Qc to be infinity. Moreover, we assume that V{r) is polynomially 
bounded at infinity, i.e, there exist constants Ci, C2 and 2 < s < 00 such 
that V{r) < Ci + C2r'^- For convenience, let inf,. V{r) = 0. 
Let i?^^(a, Q) be the ground state energy of i.e, 

E«P(a, Q) = inf e Q{Ho) H L\R'), H^lh = l} , (1.5) 

being finite for \Q\ < Qc and a > 0. Using standard methods (see e.g. [Q) 
one can show that there exists a minimizer cj)'^^ for S'^^ as long as \Q\ < flc, 
i.e., the infimum is actually a minimum. For \Q\ > Qc the functional S'^^ is 
not bounded from below. 

The purpose of this paper is a detailed study of the GP functional (p.. 31). 
One of our main results is that for any i7 > and for large enough inter- 
action strength a, no minimizer of S'^^ is an eigenfunction of the angular 
momentum L. Since S'^^ is invariant under rotation of 0, this result means 
that the rotational symmetry is broken in the ground state. This has in- 
teresting consequences on the multi-component Bose gas (or equivalently, 
Bose gas with spin). In particular, we will show that the ground state en- 
ergy (of the natural generalization of the GP functional to multi-component 
systems) depends non-trivially on the number of spin components, and the 
different components necessarily have a different wave function in the sym- 
metry breaking regime. 

The paper is organized as follows: In Section g we study stationary points 
of the GP functional, in particular minimizers of S'^^ restricted to the sub- 
space of eigenf unctions of the angular momentum with fixed eigenvalue n, 
so-called vortex states. We show that for large enough angular momentum, 
these vortex states can never be the absolute minimizer of the GP functional, 
uniformly in the coupling constant a, which will be crucial in the proof of 
symmetry breaking. In Section ^ we study the critical values of the angu- 
lar velocity Q for that an n + 1-vortex becomes energetically favorable to 
an n-vortex. These critical velocities all tend to zero as a goes to infinity, 
which will allow us to conclude that all vortex states with angular momen- 
tum smaller than a certain value cannot be the actual minimizers of S'^^ . In 
Section H we will use these results to prove symmetry breaking. Section ^ is 
devoted to the study of a GP density matrix functional, which will be useful 
in investigations on the multi-component Bose gas in Section ^ There we 
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show that in the symmetry breaking regime, the GP energy depends non- 
trivially on the number of spin components. In Section ^ we finally consider 
the special case of an harmonic potential V{r) = r^, where we derive ex- 
plicit upper and lower bounds on the critical coupling constant for symmetry 
breaking. 

2 Vortex states 

Given any stationary state of S*^^, i.e., a function with ||0||2 = 1 satisfying 

(i/o + 2a|0p - /i) = (2.1) 
for some /i G M, we define the (real) quadratic form Q{w) by the perturbation 
£f[<P + ew] - = s'Qiw) + Oie') (2.2) 

as e — >• 0, where ^^^^[0] = ^^*^^[0] — / A simple calculation shows that 

Q{w) = {w\Ho + 4a|0|2 - ^\w) + 2a^J '^w^, (2.3) 

where 3? denotes the real part. Multiplying (|2.1|) with and integrating 
shows that 

^ = £^P[^]+a J |0r. (2.4) 

DEFINITION 1 (Stability). We say that a stationary state is stable if 
and only if Q{w) > for all w G Q{Hq) fl L^(M^) that are orthogonal to 0. 



We have Q{i4>) = 0, which corresponds to a simple phase change in ( p.2|) . 
Moreover, Q{d(j)/dip) = because of rotational invariance of S*^^. Note that, 
by definition, an absolute minimizer of S'^^ is necessarily stable, since 



£^^[(0 + ew)/\\(l) + Ewh] = ^^^[0] + e'Qiw) + Oie') (2.5) 
if {(j)\w) = 0. 

We now look for special solutions to (|2.1|) of the form 



0(x) = /(r)e^"^ (2.6) 
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for some n G N, a so-called n-vortex. Here / is a real radial function. Since 
£^'^^[/e~*"'^] = £^'~^^[/e*"'^] + 2nQ we can restrict ourselves to non- negative n 
without loss of generality. At least one solution of the form ( |2.6| ) for each 
n always exists, as one easily sees by minimizing the functional S'^^ in the 
subspace of functions with L(j) = ncj). For of the form (|2.6|) there is the 
following direct sum decomposition of Q: Writing w(x) = X]m>o'^™(^) with 

w„(x) = A„(r)e*("-'")^ + 5„(r)e*("+")'^ (2.7) 

one easily sees that Q{w) = J2m>oQi'^rn)- 

Considering the stability of vortices, we restrict ourselves to the case 
Q < Qc, since for > fic all states are certainly unstable. (The case Q = Qc 
depends on the particular form of the potential V.) First of all, has only 
a chance of being stable if / has no zeros away from r = 0. More precisely, 
the following proposition holds. 

PROPOSITION 1 (Instability for /'s with zeros). Assume that either 
n = and f has some zero, or n>l and f has some zero away from r = 0. 
Then is unstable. 

Proof. For some real h we choose w{x.) = ih{r)e^'^'^ as a trial function for Q. 
We get 

Q(w) = {h\-A + — ~nn + Vir) + 2af^ - fi\h) = {h\H\h). (2.8) 

We know that Hf = 0, but because of its zeros, / cannot be the ground 
state of H, so there exists an h, orthogonal to /, with {h\H\h) < 0. □ 

Now let be the minimizer of in the subspace with angular momen- 
tum n. Then / defined in ( |2.6|) minimizes the energy functional 

£r.[f] = (/I - A + - + V{r)\f) + 2ixa / |/(r)|Vrfr (2.9) 

Jo 

under the condition 2tt J \f{r)\'^rdr = 1, with corresponding energy 

Enia) = £M] = ^^^[/e^"1 + ^n. (2.10) 

In the following, we will study Sn for all n > 0, not only for integers. Denote 

fj, = iJ, + nn, (2.11) 

which is independent of Q. The minimizer / of Sn has the following proper- 
ties. 
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LEMMA 1 (Properties of /). /(r) > for r > 0, f e C°°(R+) if V e 
C°°, and f{r) = 0(r") as r ^ 0. Moreover, f e L°°(M+), and 

II/IIL<£- (2.12) 

Proof. The regularity and strict positivity follow in a standard way from the 
variational equation for / (cf 0). Writing /(r) = r^g[r) we see that g 
minimizes the functional 



S[g] = J g{r) \^-g"{r) —g'i^) + V{r)g{r) + ar^'^g{ry j r'"^^dr 

(2.13) 

under the condition 27r J g[r)'^r'^"''^^dr = 1, from which we conclude that g is 
a bounded, strictly positive function. 

The bound (p.l2|) is proved analogously to Lemma 2.1 in 0: Let B = 
{x, 2a/(|x|)^ > fi}. We see that —A/ < on i3, i.e., / is subharmonic on 
B and therefore achieves its maximum on the boundary of B. Hence B is 
empty. □ 

We remark that all the properties of / stated in Lemma |I|, except for 
the positivity, hold for all ra-vortices and not only for minimizers. Also the 
following lemma holds true for arbitrary vortex states. 

LEMMA 2 (Properties of g). Let f{r) = r^g{r) he a stationary point of 
£n, for n e M+. Then 



gWoo <\\f Woo {clfiT^\ (2.14) 



where 



n/2 „ \ V" 



2"" (^) ^Csc (^) 
(2 - n)r(n) 

n r(n) 



for n < 1 



forn>l. (2.15) 



// / is the minimizer of En, and if V is monotone increasing, then g is a 
monotone decreasing function. 
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Proof. By a rearrangement argument one sees from ( p.l3| ) that the minimizer 
of S is monotone decreasing if V is monotone increasing. For a general n- 
vortex, g fulfills the equation 

2n -\- 1 

-g"{r) — g{r) + V{r)g{r) + 2ar'^''g{rf = fig{r). (2.16) 

Kato's inequality and the positivity of V imply that 

-\g{r)\" ^\g(r)\' < ^\g{r)\ (2.17) 

r 

in the sense of distributions. Now let Xni'f', s) be the kernel of the operator 

2n + l d , , 

+ 1 , (2.18) 



dr"^ r dr 
acting on L^(M_|_, r^""'"^(ir). It is given by 

/ X 1 fln{r)Kn{s) forr<s 
(rs)" [Kn{r)In{s) for r > s, 

where /„ and Kn denote the usual modified Bessel functions. Note that both 
In and Kn are positive, so Xn is positive. By scaling, the integral kernel of 
( |2.18| ) with +1 replaced by is t'^"'Xn{ft, st). Therefore ( |2.17| ) implies that, 
for t > and < a < 1, 

POO 

\9ir)\ < + Xn{rt,st)\g{s)\s^-+'ds 

Jo 

POO 

< (/i + t')t'1|^||L"1l/ll^ / Xnirt,st)s'-^'--'^ds. (2.20) 

Jo 

We now claim that J Xni'f', s)h{s)s'^^~^^ds is monotone decreasing in r if /i is a 
positive, monotone decreasing function. To prove this, it suffices to consider 
a step function h{s) = Q{R — s), R> 0. A simple calculation yields 

Jo I i?"+V„+i(i?)^;H forr>i?. 
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which proves the claim, since In{r)/r^ and Kn{r) / are monotone increasing 
and decreasing, respectively. Therefore the maximum on the right hand side 
of ( p.20[ ) is achieved for r = 0, and we get 

Moo<^\\9\\]^imY^J^ K^is)s-^'—ds. (2.22) 

The last integral can be evaluated explicitly, if na < 2. Choosing a = 1 
for n < 1 and a = 1/n for n > 1 and optimizing over t yields the desired 
result. □ 

Equation ( p.l4|) effectively gives a lower bound on s, the size of the vortex 
core, defined by ~ ||/||oo(^/s)" as r — 0, i.e., 

nMY"'>-(m"'>-^- (2.23) 



--Or"||/||oo; ~ Vll^lloo/ ~ fl^/^Cn 

Note that c" — > 1 as n — > 0, and c„ = 0(n^^/^) for large n. The latter fact 
will be important in the proof of the following theorem. 

THEOREM 1 (Instability for large n). For a// < < fic there exists 
an Nq < oo (independent of a!) such that all vortices with n > Nq are 
unstable. 

Proof. Let n > 1, and let Wi G i/^(]R^) be radial and normalized, with 
support in the ball of radius 1. Let X = J |wi(r)pV"(r/c„/i^/^)(i^x, T = 
{wi\ — Awi), and define w by 

w{^) = Cnfi^/^w^{cnfi^/^r), (2.24) 

with c„ given in (|2.15|) . We have, using ( |2.12|) and ( p.l4| ), 

Q{w)<nn + fl (^clT + jX-l + 2 j |w(x) |^ min{l, (r Vn)"}^^^x^ • 

(2.25) 

With M = J l^ipr^" this gives 

Q{w) <nn + jl [clT - 1 + 2M) + X. (2.26) 
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Now jj is larger than Cn^ = inf spec (—A + V) \L=n, which can be bounded 
below as follows. By assumption (|1.4| ), V{r) > — Cq + Cl'^r'^/A for some con- 
stant Cq and < < Denoting by ipn^ the eigenfunction corresponding 
to en \ we have 

ei°) > + V{r)m > m -C^ + ^ + ^rV:^) > -C^ + nn. 

(2.27) 

Now = 0{n^^) as n — i> oo, and the same holds for M. X can be bounded 
by sup^<(^2^)-i/2 V{r), which is, for fixed Q, bounded independent of n and 
a by the considerations above. Therefore Q{w) < for n large enough. □ 

For a special class of potentials, we can extend the previous result in the 
following way: 

THEOREM 2 (Instability for special Vs). Assume that for some d e 
N, d>2, 



for all r. Assume also that n > d and 



< (2.28) 



fi>nQ{l + ) . (2.29) 



d- 1 
Then is unstable. 

Proof. For 1 < d < n we choose ClS db trial function 

u;(x) = {A{r) + B{r)) e*("-°')^ + {A{r) - B{r)) e*("+'^)^. (2.30) 
Then Q{w) can be written as 



QM = 2{ ^ 



H 



A 
B 



(2.31) 



where 



H 



^0 + -nn-fx + 6a/" dil 



(2.32) 
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We now choose A{r) = f'{r)/r'^ ^ and B{r) = nf{r)/r'^. Note that A — 
B = r"-'^+i(//r")' < 0(r) as r ^ 0, so w G H'^{R^). A straightforward 
calculation using Equation yields 



/"°° f'(r)f(r) 

Q{w) = 871 dr^^-y^(-/i(rf-l) + 2a(d-l)/(r)2 + nfi 

+ {d-l)V{r)-'-V'{r)) 
firV' 



1^2(d-l)+l / ^ / ^(^ 




' " (.F^it^y ) M . (2 33) 

where we used partial integration in the last step. Estimating ci/(r)^ by 
(|2.12| ) this shows the negativity of Q{w) as long as ( p.28| ) and (p.29|) are 



satisfied. □ 

In the case of a homogeneous potential V{r) = r^, 2 < u < oo, the 
condition ( p.28| ) is fulfilled ior u = 2{d — 1), d E N, showing that in this case 
every vortex with n > d = + 1 is unstable, if (|2.29|) is fulfilled, i.e., if /2 is 



large enough. (For fixed n this means that a has to be large enough.) 

REMARK 1 (Translational stability). The calculation in the proof of 
Theorem § shows that any vortex with n > 1 is stable against translations, 
if the opposite of the assumption on V is true for (i = 1. More precisely, the 
function w(x) defined above is, for d = 1, equal to d(j)/dx. Looking at (|2.33|) 



we see that this expression is always positive for d = 1, if {rV'{r)y > 0. 
This implies that Q{d(f)/dx) > 0, and the same conclusion holds for d(f)/dy. 
Note that the condition on V is in particular fulfilled for any homogeneous 
potential V{r) = r'^. 

The choice of the test function in the proof of Thm. ^ is motivated by 
analogous considerations in [§] (see also for a treatment of the Ginzburg- 
Landau model). 
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3 The critical frequencies 

From ( p.lO|) one sees that an n + 1-vortex becomes energetically favorable to 
an n- vortex if f2 > where the critical frequency is given by 

Qnia) = En+i{a) - En{a) > 0, (3.1) 

with En{a) defined in ( |2.1CI| ). In the following we will study the properties of 
the flnS, in particular their behavior for large a. This will be important in 
the proof of symmetry breaking in the ground state of the GP functional. 

LEMMA 3 (Relation between n„'s). For n > 

^n. (3.2) 

Proof. Using fn+i, the minimizer for £n+i, as a trial function for £n+2 and 
Sn, respectively, we get 

n^+i < (2n + 3)j ^^^^l^ci^x (3.3) 



and 



> {2n + 1) / ll^+iP^d'^. (3.4) 

□ 



THEOREM 3 (Bounds on For all neNo 



O-rre 

nJa) < (2n + 1) El (a) ( 3 + 

a 



In 



(3.5) 



^n(a) > (2n + i; 



4C^ + E„+i(a) 



27re2, 

for allO<h< n^. (3.6) 



Proof. The concavity of En{a) in implies that the right and left derivatives 
of En with respect to n exist, and from the existence of a unique minimizer 
for Sn for all n we conclude that En is in fact different iable in n. Therefore 
we have 



dEn 

On 



(3.7) 



n=no 
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for some uq G {n,n + 1). Now let /(r) = r^g{r) be the minimizer of Sn- To 
obtain the upper bound, we estimate, for < a < 1, 



OK 
dn 



In 



< 2n 



< 2n 



1/2 

r-Vx ) + \\g 



x|>_R 



x|<_R 



oo^ + ||5'||ooll/lloo ^ 



for all i? > 0. Optimizing over R yields 



dn 

and by (|2l4l ), 



l/(2na+l) 



dn 

Next we choose 



<(2„ + i) (2 



1+no 2n«|| j-||2no+2 ~on\ l/(2na+l) 
^ '^n ll/lloo A'" J 



a = min < 1, 



n 



In 



4vre2||/|| 



H -1 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



and use ( |2.12| ), which yields 



dEn vre _ , , / 

< — yU max < 2?T, + 1, 77, In I — 



dn a 



cla 



27r 



Using /i < 2En{a) and c^" < e this gives, together with (|3^ ), for i7o 



27r6 r / 

^o('j^) ^ Ei{a) max |3, 1 + In ^ 



a 

2^ 



(3.12) 



(3.13) 



Now r2„ < (2n + 1)^0 by Lemma ^, which finishes the proof of the upper 
bound. 

To obtain the lower bound, we use (13.41) and 



f2 

Jn+l 



1 



> 



because of (O), for a\\ Q < Vt < Vtr. 



(3.14) 
□ 
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Note that since V{r) < Ci + C2r^ for some 2 < s < oo by assumption, 
a simple trial wave function shows that En{a) < 0(a*/*^''+^^) as a ^ oo, 
implying that f2„ behaves at most as a"^/^*^^-* In a for large a. In particular, 
lima^oo ^n(a) = for all n. 



4 Symmetry breaking in the ground state 

We now have the necessary tools to prove symmetry breaking. With the 
results of Theorems |1| and ^ the following is easily shown. 

THEOREM 4 (Symmetry breaking). For all < Q < Qc there is an 
such that a > implies that no ground state of the functional is 
an eigenfunction of the angular momentum. 

Proof. Fix < r2 < Qc- From Thm. ^ we see that there exists an 
independent of a such that all vortex states with n > Nq are unstable, 
and therefore cannot be minimizers of S'^^. By the definition of the critical 
frequencies ( |3.1| ), 

min {En{a) ~ nQ} > mm {En{a) — nQ} (4-1) 

0<ra<Afn n>NQ 



if 

0<j<Nn - j 



Q > max -TT^ — : ^ ^i{a). (4.2) 



By Thm. ^ (and the remark after the proof) this can always be fulfilled for 
a large enough, so the ground state of S'^^ cannot be an n-vortex. □ 

This shows that no minimizer of the GP functional is an eigenfunction of 
the angular momentum. We can even show more, namely that the absolute 
value of a minimizer is not a radial function. To prove this, we need the 
following general lemma. 

LEMMA 4 (Fourier series of e*^(^)). Let h : [0,27t] be a measurable 
function. Then the set of Fourier coefficients of e*^ contains either only one 
or infinitely many non-zero elements. 
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Proof. Let e'^'^^'^^ = where the hnS are the Fourier coefficients of 

e*^. Let riQ = maxjr;,, hn 7^ 0} and rii = min{n, hn ^ 0}, assuming that both 
are finite. Since 

1 = le^M'.) |2 = J2 A:„e^"^ with K = J2 K+Jh^. (4.3) 

n m 

we know that kn = 5no- But /Cno-ni = hnohn^ 7^ 0, so no = ni. □ 



COROLLARY 1 (Symmetry breaking, part 2). Let < Q < Qc and 

a > aQ, and let be a minimizer of . Then is not a radial 

function. 

Proof. Assume \(f)'^^\ is radiaL With H = Hq + 2a 10^^^ — /i we have 
HcfP^ = 0. Because H commutes with L, it has eigenfunctions hn{r)e^^'^ 
with corresponding eigenvalue 0. Therefore 



>,¥,)= ^A„/i„(r)e^"^ (4.4) 



for some 7^ € C, and the sum is finite, since the ground state energies 
of H restricted to subspaces of L = go to infinity as n 00, and there 
can only be one eigenfunction for each n. Choosing some interval / G M+ 
where |(/>^^| does not vanish, we can conclude with Lemma ^ that only one 
hn is unequal to zero in /. But since the h^s do not vanish on some open 
set, this is true on all of R+. Therefore (j)^^ has to be an eigenfunction of L, 
contradicting Thm. ^. □ 

Note that Thm. § implies in particular that the minimizer of S'^^ is 
not unique (up to a constant phase), and Corollary [l| shows that even the 
absolute value is not unique. By rotating a minimizer one obtains again 
a minimizer, which is, at least except for exceptional values of the rotation 
angle, different from the original one. 

Numerical investigations |TD|, W^ indicate that the symmetry breaking 



results from a splitting of an n-vortex into several vortices with winding 
number 1. I.e., one expects that (f)'^^ has d distinct zeros of degree 1, where 
d = deg{0*~^^/ 10*^^1} for large enough r. This property was proved for large a 
for the minimizer of models similar to the CP functional |T^, [TB|, . There 
is also experimental evidence for vortex splitting M, H. 
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5 A density matrix functional 

We now introduce a new functional, which will be convenient in the following. 
Firstly, to obtain a lower bound on the critical parameter a^, and secondly, 
for studying a generalization of the GP functional to a Bose gas with several 
components, which will be done in the next section. 

Analogously to the Gross-Pitaevskii functional we define the GP density 
matrix (DM) functional as 

£°^[7] = Tr[i7o7] + a J p,{^fd^^. (5.1) 

Here 7 is a one-particle density matrix, a positive trace-class operator on 
L^(]R^), and denotes its density. The ground state energy, the infimum 
of ( |5.1| ) under the condition Tr[7] = 1, will be denoted by E^^{a,VL). It is 
clear that E^^ < E'^^ . Using the methods of @] and ||15| one can prove the 
following theorem. 

THEOREM 5 (Minimizer of For each < n < and a > 



there exists a minimizing density matrix for l\5. J\ ) under the condition Tr[7] = 
1. The density corresponding to the minimizer, denoted by p^^, is unique 
(and therefore a radial function), and each minimizer also minimizes the 
linearized functional 

S^ni] = n{Ho + 2ap^'^U (5.2) 

The uniqueness of the density results from the strict convexity of in 
p^. In general, the minimizing density matrix need not be unique. However, 
for the functional (15.11) we can show that this is indeed the case. 



THEOREM 6 (Uniqueness of 7°^). The minimizer of S^^ , denoted by 
7°^, is unique. Moreover, it has finite rank. 



Proof. Being a minimizer of ( ^.2|) , 7^^ can be decomposed as 

7^^(x,x') = J2 hkfAr)f,ir')e^^^^-'^^'\ (5.3) 

where the fk{r)e''^'^ are the ground states of Ho + 2ap^^, and the sum is finite 
because of the discreteness of the spectrum of this operator, implying finite 
rank of 7°^^. Moreover, there can be only one ground state for each angular 
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momentum, and fj{r) = r^gj{r) with gj{r) bounded and strictly positive. 
Therefore 

p^^(r)=J2r'Xjir,^l (5.4) 

j>0 

with 

Xj{r, ^) = J2 )^j-km-kir)9k{r)e'^'-"'^^. (5.5) 

A; 

Hence each Xj has to be independent of cp as r —>■ 0, which imphes, together 
with Lemma ^ that Xj^k = for j ^ k. Moreover, Ajj is determined by the 
unique density = X]j>o^^''-^ij5'i' unique. □ 

Analogously to minimizers of the GP functional, the DM density has the 
following properties. 

LEMMA 5 (Properties of p™). For r > we have > and p°^ G 
zfVe Moreover, ||p°^||oo < p°^/(2a), where zs the chemical 
potential of the DM theory, which is the ground state energy of S^.^. 

Proof. Note that p^^ = Xjjfji^r)"^ with the notation of the proof of Theo- 
rem where Xjj > 0. The first two properties follow from this decomposition 
and a bootstrap argument. Moreover, a direct computation gives 

-Ap°^ < 2p°^(p°^ - 2ap°^ - V{r)). (5.6) 

Since V{r) > this implies that 2ap^^ < fi^^ by a subharmonicity argument 
as in Lemma |l[ □ 

An important consequence of the uniqueness of the minimizer of is 
the following corollary. 

COROLLARY 2 (Non-equivalence of and S^^). Assume that the 
minimizer of S'^^ is not unique, which is in particular the case for a > an. 
Then E^^{a,n) > E^^{a,n). 

Proof. This follows immediately from the uniqueness of j^^. □ 
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Note that in the case of non- uniqueness of </)*^^, the rank of is always 
greater or equal to two, and therefore the ground state of Hq + 2ap^^ is 
degenerate. This holds in particular in the whole region a > a^, not only 
for isolated points or lines in the (a, fl) plane. For the non-rotating case, i.e. 

= 0, E'^^ and are equal for all a. This remains true, if Q is not too 
large. 

PROPOSITION 2 (Equivalence of and ^™ for small Q). Assume 
that 



for some Q < Q < Qc- Then E'^^{a,fl) = E^^{a,fl) , and the minimizer of 
has zero angular momentum. 

Note that this Proposition implies a lower bound on a^. 

Proof. Assume that ■?/' is a ground state of = Hq + 2ap°^ with angular 
momentum Lip = mip. Then 

2 

Tfl 

f/"^ = mi spec H""^ = {\iIj\\- A -Qm + — + Vir) + 2ap''^\\iP\) 

2 

> p^^ + iijl'^-nmli;), (5.8) 



implying that 



m < — < ^(V^lr^lV;). (5.9) 



Choosing some n with n<n<ncwe have < 4(C^ + V{r))/n^ by ( |OD , 
and hence 



1-2^1 (^Ir^i 



< + (5.10) 



where we have used that the ground state energy of —A — QL + is 
Q. Therefore 

fi(^kV) < ^ < 1 (5.11) 
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if 



^ 1 . ^ 



showing that any ground state of necessarily has angular momentum 

m = 0. □ 



6 The multi-component Bose gas 

We now consider the Gross-Pitaevskii theory of a rotating Bose gas with 
different components, or equivalently, a Bose gas consisting of particles with 
spin {uc — 1)/2. The natural generalization of the Gross-Pitaevskii functional 
is 

Tic p 

S^^[^^,...,^r^^ = J2{^^\Ho<P^)+a ^ /|0d'l0,f, (6.1) 

i=l l<*J<"c 

which has to be minimized under the constraint Yli=i / — 1- The cor- 
responding ground state energy will be denoted by E^^{a, Q). 

Using standard methods, one can show that for all values of ric G N, 
< Q < Qc and a > there exist minimizing functions . . . ,0^5 for 
Sn^- The proof goes analogously to the proof of Thm. |, noting that 
can be considered as the restriction of to density matrices of rank less 
than or equal to ric- However, since this set is not convex, we can in general 
not conclude that the density of a minimizer, ^"l^ I^P^P? is unique, as it 
was the case for the DM functional. 

We see that for all values of a, and ric we always have E^^ < E^^ < 
E^^. Denoting 

n°^(a,fi) =rank7°^, (6.2) 

which we showed to be finite in Thm. ^, we can distinguish the following 
cases. 

THEOREM 7 (Minimizers of the mult i- component GP functional) 

Let . . . , (pn^ be minimizers of . 
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(i) If nc> n 



DM 



then E^^{a,n) = E\ 



(a, VL) . Moreover, 



,DM 



(6.3) 



i=l 



implying that the 0^ 's can he written as 0^ 



,GP 



Z]j=i ^iji^j' where 



^DM _ \ipi){i/ji\, ipi orthogonal, and A is an ric x -matrix 

with A^A = 1. 

(ii) If ric < n^^, then E^^{a,Q) < E^^{a,Q). 

(Hi) If Tic > 2, a > d^, then E^^{a,Q) < E'^^ {a,Q) , and the minimizers 
(f)f^ are not all equal; i.e., Yl^=i \(p?^){4>?^\ has at least rank 2. 

Proof. As remarked earlier, E^^ < E^f . To prove (i), we write 7°^^ = 

X^iLi where the are orthogonal (but not necessarily normahzed). 

Using 0j = ipi for i < and (pi = for < i < ric a.s trial functions, 
we obtain E^^ > E^f . Since is unique, ( p. 3D follows, (ii) is a trivial 
consequence of the uniqueness of 7^^. 

For a > a^, we know from Corollary [l| that there are at least two different 
minimizers, (p^^^ and (f)^'^\ for £'^^ , whose absolute values are not the same. 
For ric > 2 we use as trial functions (pi = (p^^\ (pi = (p^'^^ for i > 2 to obtain 
E^J < E'^^. Therefore the minimizers for S^^ cannot be all equal, and (iii) 
is proved. □ 

An important consequence of part (iii) of the theorem above is that 
the GP ground state energy depends non-trivially on the number of spin- 
components, at least in the symmetry breaking regime a > a^. Moreover, 
there is a clear separation between different spin-components, their individual 
densities can never be all equal. 



In the case of a harmonic potential V{r) = the theorems above are more 
explicit, setting = 2 and Cq = 0. Moreover, the special case Q = Qc = 2 
is easy to treat: only for a = there is a minimizer for S'^^ (in fact there 
are infinitely many), whereas for a > there is no minimizer, and also all 
n-vortices are unstable. 



7 The special case V{r) 
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The following bound on the energies En{a) can be easily obtained, and 
will be used several times in the considerations below. 

LEMMA 6 (Upper bound on ^„(a)). 



En{a) < 2{n + 1) /l + (7.1) 
V bn{n + l) 

with bn = 27r4"(ra!)V(2n)!. 

Proof. This follows from a trial function of the form Cr'^ exp(— cr^), where 
C is a normalization constant and c is to be optimized. □ 

The estimate above implies that, for ft the chemical potential correspond- 
ing to the minimizer of Sn, 



fi-2n< 2{En - 2n) - 2 < 2 I 1 + 2Wa ^"^ ^' I . (7.2) 



This will be useful for an upper bound on /, since in the case of V{r) = 
we can improve the estimate ( |2.12| ) by 

l<^(fc-2n). (7.3) 



The proof is analogous to the one of ( |2.12| ), using in addition that n'^/r'^+r'^ > 
2n. 

We now consider the stability of n- vortices 0, i.e., solutions to ( p.l| ) of 
the form (|2.6|) . In addition to the results of Section ^ we can state another 
proposition. 

PROPOSITION 3 (Instability for small a). Assume that a < nn{2 - 
Q) . Then cj) is unstable. 



Proof. Let n > 1, and let w{'x) = a/TTtt exp(— r^/2) be the ground state of 
Hq. Then 

Q{w) = 2 - /i + 4a j w^\(j)\'^ < n{n - 2) + ^, (7.4) 
where we used Cauchy-Schwarz and the fact that fi > 2 — n{Q — 2) + 2a J 

□ 
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Moreover, we can improve Theorem ^ in the following way. 

THEOREM 8 (Instability for large n, harmonic potential). Assume 
that n > 10 and 

where dn is a monotone decreasing function of n, with dn < 1 for n > 10, 
namely 

rf. = min (4 + ^""^i^^.^.F + 2^-" (n! - T{n + 1,2)),^]. (7.6) 

Then is unstable. 

Note that since fl > 2n, (|7.5|) is in particular fulfilled if (i„ < 1 — ^7/2. 

Proof. Let n > 1, and let Wi G if^(M^) be radial and normalized. Let 

T = {wi \ — Awi) and X = {wi\r'^Wi), and define w by 

w(x) = Cn/i^/^wi(c„/i^/V), (7.7) 
with c„ given in (|2.15| ). Using ( |7.3|) and ( |2.14|) , we can estimate 

Q{w) < nn + fi(clT + ^X-l 

+2{fl-2n) j |u;(x)pmin{l, (r2/i4)"}dV (7.8) 

With = j^^^ \ and M = J^^^ jwipr^" this gives, using /2 > 2(n + 1) in 
front of X, 

Qiw) <nn + fi [clT - 1 + 2 (AT + M)) + (^-^^^i-^X - An{N + M)^ . 

(7.9) 

Now if we choose wi{r) = (2/7r)^/^ exp(— r^) the last term in (|7.9| ) is negative, 
and the computation ofT, N and M yields 

Q{w) <nn- fi{l~dn), (7.10) 
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where dn is the first part in the parenthesis in ( [7.6| ). For large n, this can be 
improved by choosing wi{r) = (35/97r)^/^[l — r^^_|_, which gives 

Q{w)<nn-fi{l-—]. (7.11) 



n ^ 

□ 

In |ll| the authors used a similar method to the one of Thm. |^ and a 
particular assumption on the form of the vortex state to obtain a dn in 
( |7.5| ) that is less than 1 if > 2. 

We conjecture that in the case of an harmonic potential an n- vortex with 
n > 2 is unstable, for all values of > and a > . However, we can prove 
this only for Q < 1. Namely, if we insert V{r) = in ( p.33| ), set d = 2, use 
the improved bound ( |7.3| ) and /i > 2 — n{Q — 2) + 2a J |0|^ this shows the 
negativity of Q{w) as long as n > 2 and 

f^<i+^(i+«y"i0r) (7.12) 

is satisfied. This implies that all vortices with n >2 and a > are unstable 
as long as < 1. 

As a consequence of the considerations above we can state an explicit 
condition on a where an ?7,-vortex is necessarily unstable, using the general 
lower bound 



/ 



which can easily be proved using elementary calculus of variations. Moreover, 
since in two dimensions / |0|^ scales as 1/ (length)^, the virial theorem implies 
for the minimizer / of Sn 

f{rYr'd'^ = ^E^{a). (7.14) 

Hence we only need an upper bound on En{a), which is given in Lemma |^, 
to obtain a condition on a for validity of ( [7.12|) . 

The critical frequencies fin (a), defined in (^.Ij), have the following prop- 
erties. 
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LEMMA 7 (Properties of critical frequencies). For all n E No we have 
^n(O) = 2 and for all a > lim„^oo ^n(fl) = 2. Moreover, 

n'(0) = i TT^^^ < (0) < 0. (7.15) 

Proof. The first assertion follows from -E'„(0) = 2(n + 1). Using the harmonic 
oscillator eigenstates Xn{r) = ^yi/nnlr'^ exp(— r^/2) as trial functions the 
second assertion is proved by 

n„<2 + a/|x„«r and n„ > 2 - a/ |x„|'. (7,16) 



noting that J \xn\^ = 0{n ^^'^) as n ^ oo. To prove ( 7.15| ) we use the 
Feynman-Hellmann principle to calculate 

K{a) = I l/.+ir- / \fn\\ (7.17) 
where is the minimizer of Sn- For a = we have = Xn? yielding 

dZH). □ 



In the special case of a harmonic potential, the results of Thm. ^ can be 
improved. We get the following bounds on the critical frequencies. 

THEOREM 9 (Bounds on harmonic potential). For all n e No 

fi„(a) < (2n + 1)^ (^1 + j (3 + [in (^)] J , (7-18) 

^ / N 2n + 1 , , 

l^n(a) > , (7.19) 

wzt/i bn given in Lemma 



Proof. We proceed as in Thm. ^, but now use the improved estimate ( [7. 3D 
to replace ( |3.12| ) by 



dn a 1 ' V 27r 



<— (/i-2n)max<^2n + l,nln^U. (7.20) 
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Inserting ( |7.2| ) and using {n + < (27r) ^ for < n < 1 and < e this 

gives for Qq 



f!,(a)<H=£(l + yfj,„ax{3.1 + ln(iL)} 



(7.21) 



and using Lemma ^ we obtain ( [7.18|) . 



For the lower bound we proceed as in (|3.14|) and note that / 



En+i{a)/2 by the virial theorem. (|7.19|) is obtained by inserting the bound 
(O) for En+i{a). □ 



Numerical investigations in |]T6| indicate that for all n > 0, fi„ is strictly 
monotone decreasing in a, and, for a > 0, Qn < ^n+i, i-e., En is convex in n. 
Note that the theorem above states that f2„ behaves at most as a~^/^ In a for 



large a, in accordance with previous considerations (see |rT| and references 
therein) . 

We can now use the results above to derive explicit upper and lower 
bounds on a^. Denote 



27re 



^ a 




In 



(7.22) 



which is a strictly monotone decreasing function of a. We know from Thm. 
^ and ( [7. 121 ) that all n- vortices are unstable for n > Nq, where 



2 ioTn<l 
^ forfi>l. 



(7.23) 



Using the bound on the critical frequencies ( |7.18| ), we see by analogous con- 
siderations as in the proof of Thm. ^ that symmetry breaking occurs if 



E a) < 



2Nn-l 



(7.24) 



From Thm. | we see that E^^{a,Q) = E^^{a,Q) iiQ< l/fi^^. More- 
over, it is easy to see that the same holds for Q < 2 — a/vr. Namely, using 
Cauchy-Schwarz and the fact that / (p^^)^ is monotone decreasing in a (be- 
cause of concavity of E^^ in a) we have, for = Hq + 2ap^^, 



inf spec H^^ \l=o< 2 + 2a / p°^x^ < 2 + 2a 



X 



a 

2 + -, 

TT 



(7.25) 
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(7.26) 



showing that, for a < 7r(2 - f2), H^^ has a unique ground state with zero 
angular momentum which necessarily also minimizes the GP functional. 

Since the minimizer of the GP functional is unique and therefore an an- 
gular momentum eigenfunction as long as = E'^^ , we obtain as a con- 
sequence a lower bound on a^, using < 2E^^{a,Q) < 2Eo{a) and the 
bound on Eola) given in ( |7.1| ). Thus we have proved the following Theorem: 

THEOREM 10 (Bounds on an). In the case of a harmonic potential, the 
critical parameter for symmetry breaking fulfills the bounds 
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